Initial decoherence in solid state qubits 
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We study decoherence due to low frequency noise in Josephson qubits. Non-Markovian classical 
noise due to switching impurities determines inhomogeneous broadening of the signal. The theory is 
extended to include effects of high-frequency quantum noise, due to impurities or to the electromag- 
netic environment. The interplay of slow noise with intrinsically non-gaussian noise sources may 
explain the rich physics observed in the spectroscopy and in the dynamics of charge based devices. 
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Considerable progress has been recently achieved in 
implementing qubits with superconducting nanocircuits. 
Coherent oscillations 00,0] and entanglement of coupled 
charge qubits^] have been observed. Limitations in the 
performances arise from noise due to material and device 
dependent sources El 0, IS El E3 ■ Noise due to individ- 
ual impurities behaving as bistable fluctuators (BF) is a 
severe source of dephasing for charge based devices. Sets 
of BFs determine 1/f noise ^3 01 1 and effects due to in- 
dividual BFs has been observed both in spectroscopy and 
in time resolved dvnamics|lrl IbH]. Observations show a 
variety of features, as the drastic reduction of the ampli- 
tude of the coherent signal [H EL Eil] or relaxation limited 
decoherence 0, strongly dependent on the particular de- 
vice and on details of the protocol^, Il5l Il6| , Theories 
BFs environments^ IE Ell E3 allow to understand 
several physical aspects, although a quantitative frame- 
work embedding the variety of phenomena is still miss- 
ing. Phenomenological models of the environment as a 
suitable set of harmonic oscillators^ 01 E3 have also 
been studied. While they are unable to describe aspects 
related to the discrete nature of noised 0, 0], 
gaussian environments may sometimes provide useful in- 
formation. 

In this work we study numerically a model of discrete 
noise which potentially explain the experimental features 
due to 1/f noise, and seek a classification of the possible 
effects on the basis of simple theoretical arguments. In 
particular we study inhomogeneous broadening due to 
slow noise and its interplay with additional noise sources, 
pointing out that the presence of BFs may pose reliability 
problems for charge based devices. 

We consider a qubit anisotropically0 coupled to clas- 
sical stochastic process The Hamiltonian is 

H = H Q -±ffl* s (I) 

where Hq — — refers to the qubit. Both the operat- 
ing point, i.e. the angle 9 between z and f2, and the split- 
ting are tunable. This also modulates sensitivity to 
noise. For weak coupling the relaxation Tj -1 = s 2 S(Q) /2 
and the dephasing rate T 2 _1 = (2T1)- 1 + T 2 _1 , T 2 _1 = 



c 2 S(0)/2 being the adiabatic term which gives secular 
broadening 18J, are tuned by c = cosr? and s = sin 6. 
Only the power spectrum of noise, S(u>) — en- 
ters therefore in weak coupling the qubit is sensitive 
only to properties of the environment at the level of 
two point correlations. This picture breaks down if 
the environment extends to low frequencies For in- 
stance Random Telegraph Noise (RTN) due to a single 
BF, = {0, wo}, switching at a rate 70 is slow if 

go = (O' — fi)/7o > f 7], where the qubit frequencies 
fl and f2' = £1[(vq/SI + c) 2 + s 2 ] 1 / 2 correspond to the two 
values of £. This model describes an incoherently switch- 
ing charged impurity close to a qubit. For go > 1 features 
of the discrete nature of the BF become apparent • 

A set of Nbf BFs (&) switching at rates ji, coupled 
with the qubit via £(t) = Si£i(*)> models 1/f noise if 
7 2 are distributed [lit with « V7- Thc BF" 1 // 

spectrum is S 1 /^ = X)i | vf 7*/ (7 2 + uj 2 ), and if 7 l G 
[lm,lM], in the same interval of frequencies is approx- 
imated by S 1 //^) ~ [(^4) N bf v s /hi(<y M />y m )] or 1 . 
Noise extends for several decades and in particular slow 
BFs (gi < 1), an environment with memory, make unsta- 
ble the calibration of the device. Hence the qubit dynam- 
ics will depend on details of the protocol. Decoherence 
due to BFs 1/f noise for various protocols has been stud- 
ied for = 0, where exact solutions are available 0. On 
the other hand the splitting is less sensitive to fluctua- 
tions at optimal working point 0, 9 = ir/2 (parameters 
gi become smaller), and part of the effects of the slow 
noise is eliminated (at lowest order T 2 _1 vanishes). 

Ideal quantum protocols assume measurements of in- 
dividual members of an ensemble of identical (meaning 
that preparation is controlled) evolutions of the qubit, 
defocusing occurring only during the time evolution. In 
practice for solid-state devices one collects several qubit 
evolutions, in an overall measurement time t m . Lack 
of control on the environment preparation determines 
defocusing of the signal, analogous to inhomogeneous 
broadening in NMR 18]. This is also true for single-shot 
measurements!^. In our case BFs active in additional 
broadening have 7* > 7* ~ min-IWlO. t~ 1 1 [l6j . 

We first study Hamiltonian by simulating the 
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FIG. 1: Simulations of an adiabatic BF — 1/ f environment 
at 8 = 7r/2. Relaxation studied via (o x ) (green line) is well 
approximated by the weak coupling theory (dots) . Dephasing 
in repeated measurement damps the oscillations (thin black 
line) . Part of the signal is recovered if the environment is re- 
calibrated (thin gray line) . Noise is produced by = 250 BFs 
per decade, with l/t m — 10 5 rad/s < 7i < "/m = 10 9 rad/s 
<Q — 10 10 rad/s. The coupling v = 0.02 SI is appropriate to 
charge devices, and corresponds to S(uj) = 1&-kAEq/uj with 
A ~ 10 -6 01- The adiabatic approximation Eq.© fully ac- 
counts for dephasing (red dot-dashed line). The Static Path 
Approximation (SPA) Eq.© (blue solid line) and the first 
correction (blue dashed line) account for the initial suppres- 
sion, and it is valid also for times t 3> 1/tm- In the inset 
Ramsey fringes with parameters appropriate to the experi- 
ment [J (thin black lines). The SPA (blue solid line), Eq.@, 
is in excellent agreement with observations [T^|. and also pre- 
dicts the correct phase shift of the Ramsey signal (blue dots, 
compared with simulations for small detuning S — 5Mhz, 
violet line), which tends to ~ 7r/4 for large times. 



Stochastic Schrodinger Equation for the qubit in a BF- 
1// environment. We generate £(i) as a sum of Nbf < 
2000 RTN processes with proper distribution of param- 
eters. In order to minimize errors in generating £(t) we 
use a "waiting time" algorithm [2l|, which also reduces 
the computational time. The qubit propagator is evalu- 
ated as the product of the propagators between succes- 
sive switches. Finally we perform the statistical average. 
We study an ensemble of time evolutions of the qubit, 
each lasting for a time t. During the overall time t m of 
the protocol the environment evolves in an uncontrolled 
way, so BFs with 7j ^ l/i m average, whereas BFs with 
7i <C l/t m are frozen. Thus for the simulation we con- 
sider > 10 5 realizations of £(t'), for < t' < t. For 
the individual BFs at t' = we choose the same initial 
£i(0) — 0, 1 if 7j < l/t, n whereas if > l/t m we take 
a distribution with < £i(0) < 1. This prescription has 
been checked against more accurate ones in Ref.|l6(. 

Results at 9 = tt/2 for an adiabatic 1// environ- 
ment, 7m "C ^, show the presence of several time scales 
(Fig. nj. Coherent oscillations of (a y ) are initially sup- 
pressed with a power law. Relaxation occurs on much 
longer time scales, given by the weak coupling result. 
The initial suppression is due to inhomogeneous broad- 
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FIG. 2: Signal decay at non optimal points, (a) The SPA 
(solid red lines) is compared with the quadratic approxima- 
tion Eq.10 (dashed blue lines) at different bias points. Eq.© 
works near 8 = and 8 = n/2. The exact result at 8 — 
is also indicated (dots). This result multiplied by cos 2 8, is 
often used also at 8 ^ for interpreting experiments ( "Saclay 
theory" , diamonds) . Out of the optimal point the SPA agrees 
with results of the simulations with BF-1// noise, for the set 
of Fig. 1, at 8 = 5tt/I2 (b) and 8 = tt/4 (c). 



cning. This is apparent if we compare with results with 
a feedback protocol simulated by resetting £(0) at the 
same value for each realization of £(£')• 

Negligible relaxation allows to treat £(i) in the adi- 
abatic approximation. Observables are then given by 
path-integrals over a weight -P [£(£)] of the stochastic pro- 
cess. We study the averaged phase shift <£>(£), defined as 



<*(*)= \V£{t') P[£(i')] e~ l ti dt ' n W» (2) 
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which gives the decay of the qubit coherences, (a y 
exppM>(i)]. Here = n[(€(t)/fl + cf 

is the instantaneous qubit splitting. Numerical evalu- 
ation of the path-integral Eq.© fully agrees with the 
simulations. Further insight is obtained by approximat- 
ing Eq. @. The Static-Path Approximation (SPA), 
£(i') = £o accounts for lack of control on the environ- 
ment preparation via a statistically distributed £o- This 
blurs of the overall signal, an effect analogous to the rigid 
lattice line breadth in NMR [H}. For a set of BFs, if 
Nf,f is large enough £o is gaussian distributed with vari- 
ance (j^ — v 2 Nbf/4: — J °° (duj I V) S (ui) , where it is in- 
tended that we consider only active BFs, 7, > 7*. The 
result, plotted in Figs lll2l accounts for the initial sup- 
pression of the signal, showing that this latter is entirely 
due to inhomogeneous broadening. Therefore the analy- 
sis of the initial suppression may give information on the 
the amplitude of 1// noise at intermediate frequencies 
l/t m <v<l/t. A good estimate of the SPA is obtained 
by a quadratic expansion of Jl(£o) * n £° ( see FigGt) 



-i $(t) = v 14-; n 

V ' 2 1 + is 2 a 2 t tt 2 
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which is accurate close to = and 8 = n/2 for cr^/fl 
small enough. The resulting suppression factor exp(3<I>) 
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turns from a exp(— ^c 2 a^t 2 ) behavior at 9 w to a power 

law, [1 + (s 2 <7 e 2 V^) 2 ]" 1/4 , at9» 7r/2. These limits re- 
produce known results for gaussian 1// environments, 
namely at 9 — the i <C 1/7M limit of the exact re- 
sult |23| and for 6 — tt/2 the short-time result of the 
diagrammatic approach of Ref. [9|- This is not surpris- 
ing since the SPA does not require knowledge about the 
dynamics of the noise sources, provided they are slow|23|. 

Eq. can be systematically approximated by sam- 
pling better in [0,i\. For the first correction .P[£(t)] 
is approximated by the joint distribution P(£t£;£oO), 
where & = £(*)■ At 9 = tt/2 for generic gaussian noise 
we find 



i$(t) = -ln 
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where ?r(t) = ^ / °° (dw/n) S(oj)(1 - e~ luJt ) is a tran- 
sition probability, depending on the stochastic process. 
For Ornstein-Uhlembeck processes it reduces to the re- 
sult of Ref. [l^. The first correction suggests that the 
SPA, in principle valid for t < 1/jm, may have a broader 
validity (See Fig2J. For 1// noise due to a set of BFs 
it is valid also for t » 1/7M, if 7Af < ^- Of course the 
adiabatic approximation is tenable if t < T\ — 2/S(fl). 

The main effect of faster BFs in the 1 // spectrum is the 
decrease of T\ . Relaxation is due only to the fast part of 
the spectrum u ~ f2 and well reproduced by the Golden 
Rule. This is not true for decoherence, for instance in 
our example (Fig. EJ T 2 « T\, as observed in NEC ex- 
periments 0. We study the interplay of fast and slow 
noise by a two-stage elimination. We first decompose 
f (*) -►£(*)+£/(*)• Here £(*) represents all BFs having 
switching rates small enough to be treated in the adia- 
batic approximation, 7, < j a d (in practice we may take 
lad ~ O/10). Fast BFs are described by £/(t) or better 
modeled by a set of quantum impurities as in Ref. 0. 
In this case £/(<) — > £/ and we have a quantum envi- 
ronment able to produce also spontaneous decay. The 
reduced density matrix of the qubit can be written as 

P (t) = Jvat) pim Pf [t\at)] 

where i s the qubit density matrix resulting 

from the elimination of the fast environment, under the 
"drive" £tt) , and can be found within the weak coupling 
theory |24j . This is very simple if we treat slow noise in 
the SPA, where £(t) = £o- We are left with averages over 
P(£o) of the entries of p/[i|£o]- For instance the decay of 
the coherences at 9 = tt/2 is given by 



■S/(o)t-i in 1 1+ [in+s f (o)~±s f (n)] a^t/n 2 



(4) 



where Sj{uj) refers to the set of fast BFs, whereas 
refers to the set slow BFs. Eq.(0} agrees very well with 
simulations (Fig[3J). 
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FIG. 3: Results of simulations with an adiabatic plus fast 
BF — 1/ f environment (same parameters of Fig. 1 except for 
jm = 10 11 rad s _1 ). Relaxation (thick solid green line) is 
given by the weak coupling result (dots). The initial sup- 
pression of the oscillation amplitude is partially removed by 
a feedback protocol (shaded curves) and is well described by 
the two-stage elimination SPA theory (solid red line), Eq.Q. 



We notice that the validity of Eq. is not limited to 
the ~ I/7 distribution of switching rates giving rise to 
1// noise. According to this description relaxation and 
inhomogeneous broadening are due to separate sets of 
BFs. Therefore no special relation is expected to hold 
between T\ and The mixed term in Eq.Q, due to 
the interplay between slow and fast BFs does not quali- 
tatively change this conclusion. Finally Eq.Q is rather 
independent from the nature of the noise sources and the 
form of the spectrum and it is applicable in other sit- 
uations, e.g. when slow impurity noise combines with 
fast electromagnetic noise. Eq.(0} becomes exact if £f 
determines white noise, a scenario recently proposed to 
fit decoherence in phase-charge qubits. 

We come now to effects of the discrete nature of noise. 
Results presented so far rely on the SPA and on the weak 
coupling theory, therefore they apply to situations where 
discrete and gaussian noise are indistinguishable. Strik- 
ing differences appear when only decoherence during time 
evolution 0, matters, or if the distribution of envi- 
ronment couplings Vi is wide , the realistic scenario for 
the solid state. We now study the interplay of 1// noise 
with RTN produced by one BF which is more strongly 
coupled with the qubit. The model for the BF is minimal: 
it is an incoherent slow fluctuator, having 70 <C 1/t <C ^ 
but vp < tt. Even if the BF is not resonant with the 
qubit |25| it strongly affects the output signal. If 90 > 1 
[7|, it determines beats in the coherent oscillations and 
split peaks in spectroscopy, which are signatures of a dis- 
crete environment. The additional BF makes bistable 
the working point of the qubit and amplifies defocusing 
due to 1// noise. Even if the device is initially optimally 
polarized, during t m the BF may switch it to a different 
working point. The line shape of the signal will show 
two peaks, split by ~ il' — fl and differently broadened 
by the 1// noise in background. The corresponding time 
evolution will show damped beats, this phenomenology 
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FIG. 4: (a) (a y ) at 6 = tt/2, fi = 10 10 Hz. The effect of weak 
adiabatic 1/f noise (light gray line) (7 e [10 5 ,10 9 ] Hz, uni- 
form w = 0.002 fi, n<j = 250) is strongly enhanced by adding a 
single slow (7/fi = 0.005) more strongly coupled (vo/Q = 0.2) 
BF (black line), which alone would give rise to beats (red 
line), (b) When the BF is present the Fourier transform of 
the signal may show a split-peak structure. Even if peaks are 
symmetric for the single BF alone (dashed line), 1/f noise 
broadens them in a different way (solid line). 

being entirely due to the non-gaussian nature of the en- 
vironment. For illustrative purposes we show results of a 
simulation at the optimal point, where 1/f noise is adi- 
abatic and weaker than the typical noise level in charge 
qubits. This picture applies to smaller devices The fact 
that even a single impurity on a 1/f background causes 
a substantial suppression of the signal poses the problem 
of reliability of charge based devices. An analytic two- 
stage elimination combining the SPA with the solutions 
for the dynamics of a qubit coupled to an impurity 0, Q 
can be developed, and will be presented elsewhere. 

Recently effects of the resonant coupling of the qubit 
with a quantum two-level system, simulating defects in 
the tunnel oxide, have been proposed to explain features 
of the dynamics of phase Josephson qubits |25| . We have 
shown that these effects are present even if the impu- 
rity behaves as a slow stochastic fluctuator. Our model 
describes a very common situation in the solid state||, 
and it is a minimal model for charge noise in charge and 
charge-phase qubits. Finally the interplay between slow 
noise and fast noise whith generic spectrum is likely to be 
important in general and can be studied with Eq. J2Jl . The 
main open questions are the accurate characterization be- 
yond phenomenology of the physics of the noise sources 
and the design of specific strategies to defeat them and 
to improve reliability of devices. 
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